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Abstract 



Using an improved "analytic ^-matrix model", we reconsider the extraction of the cr = /o(600) and /o(980) yy widths from 
^-^ 'yy — > itn scatterings data of Crystal Ball and Belle. Our main results are summarized in Tables [3] and |4] The averaged cr "direct 

width" to yy is 0.16(3) keV which confirms a previous result of Jit] and which does neither favour a large four-quark / molecule 
C~^> nor a pure qq components. The "direct width" of the /o(980) of 0.28(2) keV is much larger than the four-quark expectation but can 
CNI be compatible with a Si or a gluonium component. We also found that the rescattering part of each amplitude is relatively large 
(J indicating an important contribution of the meson loops in the determinations of the cr and /o(980) yy total widths. This is mainly 
D due to the large couplings of the cr and /o(980) to nn and/or KK, which can also be due to a light scalar gluonium with large OZI 

violating couplings but not necessary to a four-quark or molecule state. Our average results for the total (direct+rescattering) yy 

widths: V°' - 3.08(82) keV, Tj? f = 0.16(1) keV are comparable with the ones from dispersion relations and PDG values. Using the 
£S| parameters from QCD spectral sum rules, we complete our analysis by showing that the production rates of unmixed scalar gluonia 

crg(l) and G (1.5-1.6) agree with the data from J /iff, <p radiative and D s semi-leptonic decays. 

Keywords: 

yy and nn scatterings, radiative decays, light scalar mesons, gluonia and four-quark states, QCD spectral sum rules and 
low-energy theorems. 
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1. Introduction 

In previous series of papers O]-^, we have used an improved 
version of the K-matrix model originally proposed in J4J] for 
studying the hadronic and yy couplings of the cr// (600) me- 
sonLil . We found that the "direct" coupling of the cr to yy 
is more compatible with a large gluon component in its wave 
function rather than with a qq (too large yy width) or four-quark 
(too small yy width). More recently, we have extended the anal- 
ysis for studyingthe hadronic couplings of the cr//o(600) and 
/o(980) mesons |0, [3t]. We found an unexpected relatively large 
coupling of the cr to KK: \g<rK+K-\/\g<rn+ir\ = 0.37(6)0, which 
disfavours its large n-n molecule and four-quark components, 
while the large coupling of / (980) to KK: \gj K *K-\l\gfn*n-\ = 
2.59(1.34), excludes its pure (uu + dd) content. These phe- 
nomenological observations go in lines with the fact that, in the 
7 = channel, the gluon component is expected to pla y an es- 
sential role through the scalar QCD anomaly (dilaton) | TcJ IH 
which manifests through the trace of the QCD energy momen- 
tum tensor: 

1 



2 m,(l + y m )ijfiiffi + -/3(a s )G 2 



(D 
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where y m is the quark mass anomalous dimension; if/, is the 
quark field; f3{a s ) is the Gell-Mann-Low QCD /3-function and 
G 2 is the gluon field strength. 

In this paper, we pursue the test of the nature of these scalar 
mesons by studying their couplings to yy inside the energy re- 
gion below 1 .4 GeV, where new data on yy — > tttt from BELLE 



1 20] is available in addition to the old data of Crystal Ball |21 ] 
and MARK II JH. 



2. The analytic K-matrix model for nn — > rni/KK 

In so doing, we shall work with a specific analytic K-matrix 
model originally introduced by the authors in Ref. J3l, where 
one can separate the direct and rescattering yy couplings, which 
is not always feasible using dispersion relations. 
In this approach, the strong processes are described by a K- 
matrix model representing the amplitudes by a set of resonance 
poles and where the dispersion relations in the multi-channel 
case can be solved explicitly. The model can be reproduced 
by a set of Feynman diagrams, which are easily interpreted 
within the Effective Lagrangian approach, including resonance 
(bare) couplings to nn and KK and (in the original model fl|) 
4-point nn and KK interaction vertices which we have omit- 
ted for simplicity in yl and here. A subclass of bubble pion 
loop diagrams including resonance poles in the s-channel are 
resummed (unitarized Born). In a previous work fl, we have 
discussed the approach for the case of : 1 channel © "bare" 
resonance (so-called /KD 4 model), 1 channel ffi 1 "bare" reso- 
nance (K-matrix pole) and 2 channels © 2 "bare" resonances 
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and we have restricted to the S U (3) symmetric shape function. 
We have introduced a real analytic form factor shape function, 
which takes explicitly into account left-handed cut singularities 
for the strong interaction amplitude, and which allows a more 
flexible parametrisation of the nn — > nn/KK data. In our low 
energy approach, the shape function can be conveniently ap- 
proximated by0: 



c t \ s ~ Sap n v 
fp(s) = , P = n, K , 

S - cr DP 



(2) 



which multiplies the scalar meson couplings to nn/KK. In this 
form, the shape function allows for an Adler zero at s — sap 
and a pole for crop < simulating the left hand cut. 

1 channel © 1 "bare " resonance 

Let's first illustrate the method in this simple case. The unitary 
PP amplitude is then written as: 



Tpp(s) 



Gpfp(s) 



= Gpfp(s) 

sr-s- Gpfp(s) D P (s) 



(3) 



where T PP = e iSp sin 6 P /p P (s) with p P (s) = (1 -4m 2 /s) 1/2 ; 
Gp - g^pg are the bare coupling squared and : 



Im Dp = Im {-Gpfp) = -(0p P )G P f P , 



(4) 



with: (6pp)(s) = below and (8pp)(s) = pp(s) above thresh- 
old s = Amp. The "physical" couplings are defined from the 
residues, with the normalization: 



?Ip = glpp/i^n) : IV-,™ 



Expression and normalization of the amplitudes 

The amplitude y(qt, e\) + y 2 (q 2 , e 2 ) — > n(p\) + n(p 2 ) of mass m„ 
can be written in terms of the invariants Q 

A = I l A l +I 2 A 2 , 

h = Oi ■ ei) - Oi ■ <?2)(e2 • q\)liq\ ■ qi) , 

h = (ei ■ A)(e 2 ■ A)( ?1 • qi) - (e\ ■ q 2 )(e2 ■ A)(^i • A) 
-fe • qi)(ei ■ AX42 • A) 

■ qz)(ez ■ q\)(qi • A)(? 2 • A)/(^i • qi) , (8) 

with A = p\ - p2- Helicity A = and A = 2 amplitudes are 
denoted by F and G, which, in terms of partial wave amplitudes, 
read respectively: 

F = A\ — s(s/4 - ml) sin 2 6A 2 



J] (2j+i)f J0 (s)d J Q0 (e), 



nj>0 

s(s/4 - m 2 K ) sin 2 8A 2 



J] (2J+l)g J2 (s)d J 2Q (0)4/^6- 



(9) 



evenJ>2 



dL. and d^ are the usual li-functions normalized as in PDG 12311 . 
while 6 is the scattering angle between"/? and if , which can be 
expressed in terms of s and t as: 



2t + s- 2m 2 
cos 6 = I . 

y/s(s - 4ml) 

For unpolarized photons, the cross section reads: 



(10) 



The amplitude near the pole sq where Dp(so) = and Dp(s) 
£>VOo)0 - *o) is: 



Tpp(s) 



Sq-P 



2 _ Gpfp(so) 
?crP 



-£>' P Oo) 



(6) 



The real part of Dp is obtained from a dispersion relation with 
subtraction at s = and one obtains: 



fp(s) = ~ 



hp(s) -MO) 



(7) 



with: hp(s) = fp(s)L s i(s)-(o- NP /(s-o- D p))L s i(o- D p), cr NP is the 
residue of fp(s) at cr DP and: L s i(s) = [(s - 4m 2 ,) /m p ^Li(s,m p ) 
where: L\ from yj. 

Generalization to 2 channels © 2 "bare " resonances 

In QEt], we have generalized the previous case to the one of 2 
channels © 2 "bare" resonances. 

3. The yy -» im process 



do- 2 \ 4ml , , 



(id 



where the cos 9 integration should be done from to 1 for the 
neutral and from -1 to 1 for the charged cases. In the follow- 
ing analysis, we find convenient to express the charged Fc and 
neutral F^ amplitudes in terms of the 7 = and 7 = 2 isospin 
ones |f|: 



F /=0 + 



1 



V2 

F N = -^(7^-V27^ 2 ). 
corresponding to the following \nn) states: 



\nn, I = > - 
Inn, I = 2> = 



lL+n-)- l -\n^') 



I(br + + kV>) 



(12) 



(13) 



3 Here and in the following <To is negative and is opposite in sign with the 
one used in our previous works (lKj]. 



4 We use the same normalization as @|. 



We use the same convention as |4] where Fc is opposite in sign with (27 



2 



The example of one pion exchange for yy — » n + n 
The expression of the Born term amplitude due to one pion ex- 
change reads 0: 



T B = - 
2 



2(ei ■ pi)(e 2 - pi) 2(ei ■ pi)(e 2 - pi) 



(14) 



from which one can deduce the helicity amplitudes: 



Ft = - 



mh 



2 (t - m 2 n )(u - ml) ' 
a t 2 + ts - 2m 2 n t + m\ 
2 (f - ml)(u - ml) 



(15) 



Amplitudes for 1 channel © 1 resonance below 0.7 GeV 
- The isospin 1=0 channel: starting from the S wave amplitude 
in Eq. (O, we derive the amplitude T y !) for the electromagnetic 
process for isospin / = as: 



T (0) 

y 



r 



3" 

Dp 



+ a- 



P_ 

Dp 



(16) 



Here the contribution from the Born term of yy — > n + n~ is 
given by f B = 2L\ as defined in [4], a real analytic function 
in the s plane with left cut s < 0. The function fp represents 
nn rescattering; it is regular for s < Am 2 n but has a right cut for 
s > 4ml with: 



lmf*(s + ie) = (0pf P fl)(s), 



(17) 



which vanishes at s = 0. With this definition the Watson theo- 
rem is fulfilled, i.e. the phase of 7"^ is the same as the one of 
the elastic amplitude D p l in Eq. (0. The real part is derived 
from a dispersion relation with subtraction at s = for fp to 
satisfy the Thomson limit, and has a representation similar to 
the one in Eq. (0, but by replacing (fp,hp) by {fp,hp). The 
function hp is defined as hp below Eq. (0 but with L\ replaced 
by -L 2 everywhere. It vanishes at s = cr DP and fp(s) is regular 
at this point. Finally, the polynomial P reflects the ambiguity 
from the dispersion relations and is set here to V — sF y V2. 
It represents the direct coupling of the resonance to yy. The 
residues at the pole sq of the rescattering and direct contribu- 
tions to T^ ' in Eq. ( TTBT l. respectively, are obtained as: 



6cry 6<r?r 



2 GfB(s ) 
—a- 



s Fy V2 

3~-D' p (s y °^° u " ~-D'p(s o y 
from which one can deduce the branching ratio: 



icry&o-n 



1 \8<Jy\ 



2a 2 



so 



- 1 cr—*7m IP( SO) I 1 gcrn 1 lp( s 0) I 1 Gf P (sq)\ 



(18) 



(19) 



- The isospin 1=2 channel: similarly, we parametrize the 1 = 2 

(2) 

S -wave amplitude T {) by introducing the shape function f 2 : 



-(2) 



l-A/ 2 (s)' 



Ms) 



SA2 



(s - cr D1 )( 1 S' - 0- D2 ) 



(20) 



Here and in the following, we use the same normalization as in [4] and we 
use the gauge conditions: = e;av = : £, j = 1, 2. 



and obtain: 



T (2) 

y 



a 



Afj(s) 
I-A/2WJ 



(21) 



where f 2 B = and: Im/ 2 ( S ) = (e P )f 2 (s), lmf 2 B (s) = 
(9p f 2 f 2 B ){s). These amplitudes are again both subtracted at 
s = as in case of / = and one finds in analogy: 



Ms) 



h 2 (s) - MO) 



(22) 



where: h 2 (s) = f 2 (s)L s t(s) - (<r N i/(s - L sl (o- m ) 
-(o~N2/(s - o- D2 ))L sl (o- D2 ); o- m ,o- N2 are the residues of f 2 (s) 
at cr D1 , <xd2 and f 2 B {s) is defined as f 2 (s) in Eq. (|22"1 i but with L\ 
replaced by -L 2 . The cross sections for the nn and yy scatter- 
ing processes are obtained from the previous expressions of the 
amplitudes. 

Results of the analysis 

In doing the analysis for the case of elastic nn — ♦ nn scattering 
and yy — > nn below 700 MeV where the Crystal Ball 112 1J] and 
MARK II [22] data have been used, we have obtained in yj] the 
results summarized in Table |4] 



4. Extension of the yy -* mz analysis below 1.09 GeV 

In this paper, we extend the previous analysis by including nn 
and KK loops and work in the region just above the KK thresh- 
old (the minimal x 2 /ndf of our fit is obtained for = 1-09 
GeV), where the cr(600), /o(980) contributions are dominant. 

S-waves masses and hadronic couplings 
In this case with 2 resonances © 2 channels, the hadronic 
masses and couplings of the <x(600) and /o(980) and the cor- 
responding values of the "bare parameters" of the model are 
given in Tables 2 and 3 of 13[] which we shall use in the ex- 
traction of their yy couplings. The values of the masses are (in 
units of MeV): 

Mq- = 452(12) — i260(15), M/ = 981(34) - il8(ll), (23) 

and the ratios of the hadronic couplings are: 

\gcrK+K- 



0.37(6) , !^E1£_1 = 2.59(1.34) . 

\gf7fX-\ 



(24) 



D-wave mass and hadronic couplings 

In the energy region where we shall work below 1 .09 GeV, the 
D-wave contribution can be also important. However, an ac- 
curate parametrization of the D-wave contribution is not avail- 
able. Assuming that it is dominated by the /2(1270), we extract 
its complex pole position and hadronic couplings and the cor- 
responding "bare parameters" of the model from the fit of the 
7 = 0, J = 2 phase shift measured in rf24ll and I25I1 . In so do- 
ing, we parametrize, as in 0], the 1 = D-wave mi — > nn/KK 
scattering amplitudes: 



tJ=2 
T„ K =gng K 



■ 4mi 



so 



' s — 4m 



so 



2 \ 



J=2 



(25) 



Hyams 
Froggatt 





The resulting values of the bare parameters are given in Table 
[1] from which we derive the pole position in the 2nd sheet and 
the residues of / 2 ( 1270)0. Compared with the PDG data0: 



M 



h 



1.27 GeV and T 



156.9!^ MeV 



(30) 



one can notice that the pole position and the nn width are well 
reproduced. 

Table 1 : Values in GeV* (d = 1 , 2) of the bare and physical parameters of the D-wave 
s r = 2.06 so = -1.25 gn = 0.96 g K = 1.38 

M f2 = 1.27 - i0.07 g h „ = 0.48 - i0.03 g hKK = 0.28 - i0.06 



The vector meson contributions 

In the energy-region where we shall work, exchange of vector 
mesons V = p, to, K* + , K*° in the t-channel can become im- 
portant. We introduce their couplings to yy via the effective 
interaction: 



Lvny^ e -^ Vpa Yj hv7lV ^ F P- 



(3D 



Figure 1 : a) I = D wave nn phase; b) inelasticity. Experimental data are taken 
from 



where V^ v = d„V v - d v V M (resp. F^ v = d^Fy - dyF^) are the 
vector (resp. electromagnetic) field strengths; n is the pion field. 
The coupling hy is normalized as: 



where so is a pole on the negative real axis to simulate left hand 
singularity, and: 



Dj=i(s) = s R - s- ^ g 2 P C(s, nip) 

P=n,K 

The function C(s, nip) satisfies: 

( s- 4ml 

ImC(s, mp) = Pp(s) 



{ s- So ) 



(26) 



(27) 



with pp(s) = ^jl - 4mp/s the phase-space function. The real 
part of C(s, mp) can be obtained from dispersion relation: 



C(s,mp) = — 

7T 



( s - 4ml\ 2 ( s - Ami 



\ s - s ) 
s(2s - s) 



(s - s ) 2 



A - 



-Li(s) + l 



s- s 



(28) 



where the last poles are adjusted to cancel the poles at s = so 
and0 



B = 



(s -4m 2 p \ 2 

so 

(S Q-Am 2 ^ 2 
so 



so - Amp . 



Li(s )+ 1 



s + mp _ 



Li(sq) ■ 



so + Am 2 1 



s - Amp 2 



(29) 



7 Note the extra factor 1/2 in B compared to the one in (J. 



R\' ->ny) = ^h\M\ 



1 - 



.2 \ 



My J 



(32) 



Using the standard vector form of the vector propagator, the 
Born contribution to the yy — » ti^tT amplitude due to the vector 
meson exchange is: 

t b _ ah v 
Iv "16 

(ei ■ ei)[t(s -u) + m A n ] - 2s(e t ■ pi)(e 2 ; Pi) 
t-M\ 

(ei ■ e 2 )[u(s -t) + m\\ - 2s(ei ■ pi)(e 2 ■ pi) 



u - M 2 

from which we deduce the helicity amplitudes: 
ts 



(33) 



F D 



16 



+ 



t-M\ u-M\ 



ah 2 , \ t 2 + t s - 2m 2 J + mi 



t-M\ 



u 2 + us — Imlu + mi 



u - My 



(34) 



The results agree with the ones in the different literature (see 
e.g. liil). 



"There is also pole (1.244 - i0.095) GeV in the 3rd sheet. 
'Notice that the data of the inelasticity are not quite good which induces a 
relatively bad x 1 Indf = 132.6/95 = 1.4. 



The B(l + ) axial-vector contributions 

Their contributions can be introduced via the lowest order ef- 
fective coupling ll26ll : 



L Bny = 3eh B F IJV Tr(B^Q,d v 7T}) , 



with: 



( *?_ Aid 170) 

y/2 V3 



fc+(1235) 



M 1235) + 



(35) 



(36) 



and hs is normalized as in Eq. ( l32l l. One can deduce the ampli- 
tude: 



a hi 



(ej ■ e£)(t - ml) 2 - 2s(e\ ■ pi)(e 2 ■ pi) 



t-M\ 



16 

(ei • e£){u - m 2 .) 2 - 2s(e x ■ p 1 )(e 2 ■ pi) 



u - M 2 B 



, (37) 



Then, the helicity amplitudes are: 
ah 2 , r 



16 



t - M 2 u-M 2 

b, b, 



ahl \t 2 +ts-2m 2 J + mt 



16 



t-M 2 h 

bi 



u 2 + us - 2m 2 u + ml 



u-M 2 



(38) 



The Born and unitarized S-wave amplitudes 
The Born and unitarized terms can be calculated unambigu- 
ously using the effective lagrangians. Taking into account the 
t-channel exchange of pion, vector and axial-vector mesons dis- 
cussed in the previous section and shown in Fig. [2] the Born and 
unitarized parts of the amplitude given in Eq. ( fT6b generalize to 
(normalized to a): 



1 n 
1 K 



= K 



fn + fn T nn \ , ( f%T K „ 

f?T xK ) + bK \f° + f°T KK 



+ 2 b v h 

V=p,oj 
V=K**,K*° 



2 I fy + fy *nn 
A f?f,K 



s 



byhyy sua 
BG 



r 

b A h 2 A l J * 



A~a\,b\,h\ 



fy Tkti 
+ fyTKK 

f%T nK 



(41) 



where the values of hy^ are in Eq. (l40b . b njJi ... are Clebsh- 
Gordan coefficients for projecting n,p, ... exchanges on the / = 
s-channel amplitudes: 



b„ = V2/ a/3, b p = a/3/ a/2, £ w 
^ = 1/V2, b K * + =b K ,o = l/^, 
b ai = V2/ V3, = V3/ V5, 6*, 

The partial S -wave Born amplitudes read : 



1/V6, 

= 1/V6. (42) 



_77ie ai(l ++ ) axial-vector contribution 

We describe the ai(l ++ ) in the same way as the bi(l + ~) me- 
sonF°[ where h\ n is simply replaced by h ai , which can either be 
determined from the a 1 — > ny width or from the ChPT coupling 
constants [|27l 12811 : 



hl = W±±M = 0.656 GeV 2 , 
where = 92.4 MeV is the pion decay constant. 



(39) 



The size of the different radiative couplings 

These can extracted from the data and using S U(3) symmetry 
relations and read in units of GeV 1 : 

h p = 0.82, h u = 239, ^.+ =0.83, h K ,o = l.21, 
h ai = 0.81, h bl =0.65, h hl =3h bl . (40) 



5. K-matrix model analysis of yy -* below 1.09 GeV 



10 A tensor formulation of the axial-vector meson has been proposed in l28ll 
where the form of the propagator differs from the standard one. If we use this 
propagator, we reproduce the expression of the amplitude given in 1 27] where 
an extra contact term s(t\ .tj) is added in Eq. <37t . We shall see in our analysis 
that the presence of this term would decrease the strength of the direct coupling 
of the scalar resonance but the total = direct+rescattering contribution remains 
almost unchanged. This term might be absorbed by some other counter terms 
of the complete effective lagrangian. 



2Li : P = n, K ; f» 



C = sf-i + zd; /* c = -[-i + l 2 ], 



(43) 



where L\(s,m 2 p ) and L 2 (s,m n ,M VA ) are functions analytic in 
the left hand cut plane whose expressions are given in Appendix 
C of Q]. The reduced amplitudes are defined from the ampli- 
tudes in id H as: 



T 



G_ 



7T7T 2 rj O rt 2 rj 

~PT ~ 8na F aa + ^gnagnbPab + g n b^bb 
J na 



- gnagKaPaa + {gnagKb + gKagnb)Pab + gnbgKbPbb 

= T nK , (44) 



where f m is the shape function assumed to be the same for n 
and K. fj (P = n, K) is defined in Eq. (O like f P in Eq. © 
expressed in terms of the functions hp(s) - hp(Q) but with L\ 
replaced by — L? everywhere. For vector meson exchanges, the 
triangle loop function is generalized to: 



fl{s) = ^\h B v {s)-h*m 
on 



(45) 



where 



hy(s) 



fp(s) 



2M 2 v L 2 (s) ■ 



2s(s - Am 2 p ) „ 



Li(j) 



+s 



2m 2 p - Ml 



M\ 



M 2 , 



M 2 - m 



In 



M\ 




Figure 2: Born and Unitarized amplitudes Tp\ P = n, K; V = p, <±>; A 
b\, hi, a\. 



o~m 



cr D 



2M 2 v L 2 (cr D ) 



2o- D (cr D - 4m 2 p ) , 



+cr D 



2mp - My 



L\(cr D ) 



Ml 



^ i M 2 V 



My- 



In- 



p > 



(46) 



and where the functions L\(s) and L 2 {s) have been defined in 
the appendix C of fl. 

For the axial-vector mesons, h B A (s) is opposite in sign with the 
one for vector mesons given above. 

The Born and unitarized D-wave amplitudes 
Here we shall discuss two models, namely the one originally 
discussed in and a new model with a shape function inspired 
from the S -wave channel and which is expected to have a much 
better analytic property (absence of a pole in the left cut). We 
shall assume that the D-wave contribution is dominated by the 
helicity two 1 = amplitude which will be verified (a posteri- 
ori) from a fit analysis. 

- The original model has been discussed in details in 0]. For 
one pion exchange, the isospin / = yy — > nn amplitude is: 



Am 



So 



where f B is the Born term with helicity 2: 



J IX Q 



2ml 



8ml 



s — Ami s ~ 
and the unitary triangle one loop function satisfies: 



s — Am 

Im// = 0(s - Am^ppfp 3 p - 

s- s 



(47) 



(48) 



(49) 



The real part of ff can be derived from dispersion relation with 
subtraction at s = 0: 



f 



s - Ami i S — Ami \ ~ 



s- s 



3ml 



s - Ami ~i 
■4 -L? 

s- s Q 



where 

B 2 



+ -^+B 2 
s - so 



sq - Ami ~ 
2 — L i( s o) 



(50) 



m 

m|/13 
so 



2 -^ri +4Li(io) 



(51) 



- The new model is introduced to avoid the left hand pole at 
s = so °f the old model. In this case the amplitude reads: 



f n B +g 2 JX m (s-Aml)D[ 



.1=2 



(52) 



where the 2nd term is the unitarized amplitude. f B \ ne w is the 
triangle loop function due to one pion exchange: 



6{s - Am P )p P ff- 



■ Ami 



(53) 



(s - so) 2 ■ 

The real part of f B \„ ew can be derived from dispersion relation 
with subtraction at s — 0: 



fn \new 



Ami 



(s - so) 2 
s - Ami ~ 



Ami 



(s - Sq) 



3m 2 K 
A 2 ml 
(s - So) 2 



B' 



s - s 



3 si si -so) 



(54) 



with: 



B' = lim — 

j->j ds 



(s-Am*) 2- 



3ml 



Ami ~ -, \ ~ -, 

-)L\ - A(s - Aml\l} 



(55) 



We compare these two models in Fig. [3j where we can notice 
that the two models lead (almost) to the same amplitudes. From 
this figure, it is interesting to notice that there is a strong cancel- 
lation between the Born term and the real part of the unitarized 
amplitude around the f 2 {\21Q) pole, while the imaginary part 
of the amplitude is relatively small. This feature demonstrates 
that the yy total width of the f 2 is dominated by its direct cou- 
pling as expected. In the following, we shall use the old model 
for our fitting procedure due only to a chronological procedure 
of our analysis. 

The direct resonance couplings 

On the contrary the direct couplings of the resonances are 
model dependent where the polynomial f reflects the ambigu- 
ity from the dispersion relations^- We parametrize this con- 
tribution by introducing the effective photon-photon-resonance 
couplings for the S -waves J4J]: 

V2s[ - - 

(fo-y + sf^Ta-x+iffyy + sf^Tfa 



Tt=a 



P 



Dp 



(56) 



1 1 In 12911 these polynomial (subtraction constants) are related to the pion po- 
larizabilities. 




Figure 3: D-wave amplitudes in the two different models. Born term (contin- 
uous - red line). Original model |4j]: real part (dashed line - green); imaginary 
part (dotted line - blue ). New model: the same as for the original model but 
with thick lines. 





Figure 4: Direct couplings of the resonances to yy 



where the reduced amplitudes are: 

fmr = [gna(SRh - s) - gKbignagKh ~ gKagnb)fK] , 
fnf = [gnb{SRa - s) - gKaignbgKa ~ gKbgrndU] ■ 



6. Fitting yy -* irir just above the KK threshold 

Fitting procedure 

In so doing, we fix the value of the A = 2 components of the 
/2(1270) using the fa yy width given by PDG and neglect the 
A — one. We use as inputs the Set 2 and Set 3 hadronic 
parameters obtained in [3] using the largest range of hadronic 
data. Then, we perform a fit of the four direct couplings from 
the yy — > 7t°jt° total cross-section up to yfs ~ 1 GeV^J . We 
move around 1 Ge V and looks for the minimum x 2 1 ndf 
for the total cross-section, which is obtained at -\fs= 1.09 GeV, 
where x 2 l nd f = 39.5/41 = 0.96 for e.g. Set 3 of the hadronic 
parameters, which we show in Fig. [5^). In Fig. [SJ)), we com- 
pare the yy — > ttV differential cross-section at t/s = 0.97 GeV 
for sum of all partial waves with the sum of the S+D waves and 
with the Belle data. The results of the fit are given in Tables [2] 



Belle 

Crystal Ball 
I+V+A 
D 

- E+V+A+D 
All partial waves 




(57) 



s (GeV) 



The functions fo(s), fp(s)(P = n, K), Dp(s) and bare parame- 
ters o- D , s A , s Ra , g m , g Ka , s Rb , gnb, gKb are defined in CO, H] and 
given in Tables 2 and 3 of ref. ft^. Similarly, we have for the 
D-waves 0]: 



jD _ 



a 

V5 



* 2 /Tr° + 'St 



hi 



with the normalization: 

rf^yy = UaSl; 2 fM} 2 , 



(58) 



(59) 



| 35. 
in 

8 30 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
|C0S8| 



if one assumes that the A — helicity contribution is negligi- 
bleB Using the PDG value (2.6 ± 0.24) keV Q for the yy 
width, we deduce: 



\fiy 2 \ ~ 0-136 GeV~' 



(60) 



Figure 5: a) Fit of the yy -» A total cross-section < 1.09 GeV and | cos9| < 0.8) 
for the Set 3 of hadronic parameters where ^jndf = 39.5/41 = 0.96 is minimal. 
The data come from Crystal Ball (blue full circle) and from Belle (black full square); 
dotted blue (S-channel contribution); dashed green (D-wave contribution); continuous red 
(total contribution); dotted-dashed salmon (sum of all partial waves); b) Comparison of 
the effects of the sum of all partial waves (dashed-dotted salmon) with sum of S+D waves 
(continuous red) and the Belle data on the yy — > n Q nP differential cross-section at = 
0.97 GeV. 



12 If we let free the two couplings of the A = and 2 components of the 
fl(\210) in the fit, we find that the A = coupling is negligible confirming our 
assumption. We also notice that at the /2-pole, there is a strong cancellation 
between the Born and rescattering contributions (see also (4]) which justifies 
the identification of the fi — > yy total width given by PDG | 23jl to the "direct" 
width. 



and|3] 



13 We choose not to fit much above 1 GeV in order to avoid the possible 
contribution of an eventual /d(1370) and to minimize the contribution of the 
D-wave. 
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Table 2: Fitted values of the bare couplings of a and /q(9S0) in units of GeV 1 for 
(fsy) an d GeV" 2 for (fgj) using Set 2 and Set 3 of the hadronic parameters from Qj. 
The yy — * nn total cross-section has been fitted until ^5=1.09 GeV (minimum value of 
X 2 ftidf) and the differential cross-section until = 0.97 GeV. 





fay 


fay 


flay 


f 

J fay 


x'/ndf 


Set 2 


2.68 


-2.70 


0.85 


-1.18 


51.4/41 = 1.25 


Set 3 


3.03 


-3.07 


1.03 


-1.43 


39.5/41=0.96 



Comments on the results 

Our results are summarized in Tables [3] and[4] 
- The direct part of the yy widths: one can notice that the K 
exchange and /T-loop tend to decrease the direct width of the cr 
meson which is compensated by the V+A contributions, such 
that at the end, the final result is compatible with the one 0.13 
keV obtained in the case: 1 resonance © 1 channel obtained 
below 0.7 GeV JH. We consider as a final result the average 
from Set 2, Set 3 and from the one resonance © one channel 
analysis below 0.7 GeV of IjJ] for the cr and the average from 
Set 2 and Set 3 for the /o(980) (see TableHl: 



-<dir 



■■ 0.16(4) keV, Tf * 0.28(1) keV 



(61) 



For a consistent comparison with some other theoretical esti- 
mates (QSSR,...) obtained in the real axis, we translate these 
widths from the residues to the one evaluated at the on-shell o~ 
mass or Breit-Wigner mass defined in yj, [30J] when the ampli- 
tude is purely imaginary at the phase 90°: 



Re D(Ml) = 
In this way, we obtain: 

1.2(3) keV 



M°J a 0.92 GeV. 



rdir I 
a- lon-shel! 



\on— shell 



-<dir 
-fo ' 



(62) 



(63) 



which are similar with the results obtained by using a Breit- 
Wigner parametrization of the data |l|]. 
- The rescattering part of the yy widths are (in units of keV): 



1.89(81) keV . 



Y'resc 
l fo 



0.85(5) keV . 



(64) 



where we take the average from Set 2, Set 3 and from the one 
resonance © one channel analysis below 0.7 GeV of (Q]] for 
the cr and the average from Set 2 and Set 3 for the /o(980) 
(see Table |4j. One can notice that in both cases, the rescatter- 
ings are relatively large indicating the important role of meson 
loop contributions in the yy widths of the scalar mesons. The 
large couplings of scalar to meson loops are often interpreted in 
the current literature as being related to their four-quark or/and 
molecule structure. However, these large couplings to nn and 
KK are also expected if the cr, f have large gluon component 
and violate OZI rule in their hadronic decays Hoi [l2 , 131. In 
the case of the cr meson, one can notice the large effect due to 
vector mesons which is partly compensated by the one of the 
axial-vector mesons. Compared with the result 2.7 keV from 
the analysis below 0.7 GeV with a pion loop [1] , one can no- 
tice that F r J sc has been affected by the presence of the /o(980) 
when doing the fitting procedure. In the case of the /o(980), the 



Table 3: The same fitting procedure as in Tablef2]but for the physical couplings 
(in units of a X 10 -3 GeV) and for the two photon widths (in units of keV) of 
the tr and /o(980). Each contributions of different mesons exchanged in the 
t-channel and in the loops are shown explicitly. D corresponds to the D-wave 
contribution. 



Fit 



I = n+K E+V 



E+A X+V+A X+D S+V+A+D 



Set 2 

gf -2+49i 19+15i 26+19i 10+22 i 20+27i 8+24i 9+37i 

g" sc 65+99i 66+103i -23+185i 131+28i 42+109i 66+103i 42+109i 

g<°< 63+148i 85+118i 3+204i 141+50i 62+136i 74+127i 51 + 146i 



Yf 0.21 

r™ c i.6o 

2.90 



0.06 
1.71 

2.36 



dir 

8 fa 
resc 

g fo 

gf -113+15i 91+1 li 



-125+3i 
12+12i 



-4-5i 
95+16i 



ydir 

fo 
fa 

ylol 



0.82 
0.02 
0.68 



0.002 
0.48 
0.43 



0.11 
3.97 
4.66 

-109-lSi 
47+16i 
-62-2i 

0.63 
0.13 
0.20 



0.07 
2.05 
2.51 



0.13 
1.53 
2.50 



0.07 
1.71 
2.43 



0.16 
1.53 
2.67 



26-6i -75-14i 
178+14i 130+14i 
204+8i 55 



-5-4i -73-12i 
95+16i 131 + 14i 
90+12i 52+2i 



0.04 
1.67 
2.19 



0.31 
0.90 
0.16 



0.002 
0.48 
0.42 



0.29 
0.90 
0.17 



Set 3 

gf 18+29i 12+28i 27+16i 0.2+42i 17+33i 0.3+36i 5+42i 

g£ sc 58+96i 58+100i -41 + 19H 132+20i 33+109i 58+100i 33+109i 

g'°' 76+125i 70+128i -14+207i 132+62i 50+142i 58+136i 37+152i 



Y'resc 
yioi 



dir 

g fo 
„resc 

8 fa 

lot 

s fa 
L fo 

Y>'esc 

fo 
V/o 



0.13 


0.10 


0.11 


0.19 


0.15 


0.14 


0.20 


1.36 


1.44 


4.22 


1.94 


1.43 


1.44 


1.43 


2.34 


2.35 


4.70 


2.32 


2.49 


2.43 


2.67 


99+42i 


-10-8i 


-95-14i 


22-8i 


-65-1 li 


-13-7i 


-71-10i 


12+8i 


96+6i 


51+7i 


170+3i 


124+4i 


96+6i 


124+ 15i 


87+50i 


86-2i 


-44-7i 


192-5i 


59-7i 


83-i 


54+5i 


0.61 


0.01 


0.48 


0.03 


0.23 


0.01 


0.27 


0.01 


0.49 


0.14 


1.51 


0.81 


0.49 


0.81 


0.53 


0.39 


0.10 


1.93 


0.19 


0.36 


0.15 



effect of the K and of V + A are quite important 

- The total = direct + rescattering part of the yy widths are (in 

units of keV): 



T'°' = 3.08(82) keV , = 0.16(1) keV . 



7. Comparison with some other results 



(65) 



Dispersion relations 

A comparison of the total width in Eq. d65l l with the results ob- 
tained using dispersion relations [H, 27, 3j], [32|] is shown in 
Table |4] The results can only be compared for the total width, 
as the authors have not performed the (model-dependent) sep- 
aration of the direct and rescattering processes done in this ap- 
proach. Our results for the cr are in better agreement with the 



14 Some calculations based on the n and K loops can also be found in [33 
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Table 4: Summary of the fitted values of the yy decay width(in unit of keV using Set 2 
and Set 3 of hadronic parameters from [ 3] and comparison with some other determinations 
and the PDG08 value. The total cross-section has been fitted until yfs = 1.09 GeV (min- 
imum value of x~ I ndf) and the differential cross-section unitl yfs = 0.97 GeV. The best 
X 1 (ndf = 0.96 is obtained from Set 3. 





Set 2 


Set 3 


m 


[27] 


[31] 


[31] 


[8J 


[32] 


[34] 


PDG [23] 




1.09 


1.09 


0.8 


0.8 


1.44 


1.44 


1.4 




1.44 




D/2 


0.0 


0.0 






0.13 


0.26 


0.15 








ydir 
1 cr 


0.16 


0.20 


0.13 












0.010 




yresc 
1 cr 


1.53 


1.43 


2.70 
















ytot 
1 cr 


2.67 


2.67 


3.90 


1.7 


3.1 


2.4 


2.1 


1.2 






ydir 

fo 

yresc 

Vo 


0.29 


0.27 














0.015 




0.90 


0.81 


















ytot 

h 


0.17 


0.15 






0.42 


0.10 


0.13 






0.29 ± 0.08 



ones from II Ul8l 13111 . The one of the fo is in between the results 
of |H,[H and of the PDG value Q. However, the PDG value 
seems to be more consistent with our value of the direct width 
which might be identified with the one using a Breit-Wigner 
parametrization of the data. 



_Model ofh34l 

This approach of 10,(34, 35] presents some similarities with 
ours where the direct and rescattering processes can be also 
separated. However, we differ in the treatment of the hadronic 
nn — > nn/KK process as an (a priori ad hoc) background phase 
has been explicitly introduced which multiplies the resonance 
contributions to the amplitude in this approach. Besides this 
difference, the authors also have a different (philosophical) view 
on the problem, where they claim that considering the residues 
of the cr pole for the hadronic amplitude does not give a clear 
understanding of the nature of the <x, while, in our paper, we, on 
the contrary, claim that the residues in the complex plane can 
clarify its nature. For a better understanding of the quantitative 
difference between the two approaches, we look in details into 
the analysis of J3: 

- We note that the authors use the bare couplings for predicting 
the yy width of the resonances while we transform the bare cou- 
plings (which are real numbers) to the residues in the complex 
plane by analytically continuing these results to the 2nd sheet. 
In fact, we expect that working with the bare couplings for the 
wide complex cr meson cannot (a priori) be a good approxima- 
tion. 

- We also note that the authors only include the S and D-waves 
and neglect the contributions due to the vector and axial-vector 
mesons (V+A) which is a bad approximation for the direct 
width predictions (compare the last two columns of Table |3]l 
though the numerical fits of the total cros-section are similar in 
the two cases (see Fig. [5j because of the dominant contribution 
of the rescattering amplitude which remains almost constant. 

- From Table [3] one can indeed deduce within the approxima- 
tion of [34] by retaining the S+D wave contributions that the 
value of the / (980) direct width is about (0.002-0.01) keV for 
Sets 2 and 3 which agrees with the tiny value obtained in 13411 . 

- For the cr, a comparison of our result with the one of 13411 is 



less direct due to the large width of the cr and to the important 
effect of the background phase which multiplies the resonance 
hadronic amplitude. Indeed, a factorization of this phase can 
allow a 90° phase at a low on-shell mass of about 500 MeV, 
while it is about 1 GeV, similar to the one obtained from a Breit- 
Wigner paramterization, in our approach [Eq. d62ll. If one uses 
the previous inputs (bare couplings + S and D-waves), one can 
see from the value of the bare coupling obtained in yj] that the 
direct width of the cr would be about 0.02 keV for a cr mass of 
0.42 GeV which is similar to the result of [34]. The results of 
our tests agree with the ones of ll34ll but these results might not 
be realistic due the drawbacks which we have mentioned above. 



■ Mark II & Belle 
E+V+A 
D 

— E+V+A+D 

All partial waves 




0.6 0.8 

s"*(GeV) 



Belle 
- E+V+A+D 
All partial waves 




|cose| 



Figure 6: a) Predictions of the yy —* n^n" total cross-section ( ^[s < 1.09 GeV and 
|cos#| < 0.6) using the fitted parameters from yy — > n°x Q . The data come from Mark II 
and Belle: dotted blue (S-channel contribution including direct couplings); dashed green 
(D-wave contribution); continuous red (S+D partial waves); dot-dashed salmon (sum of all 
partial waves); b) Comparison of the effects of the sum of all partial waves (dashed-dotted 
salmon) with sum of S+D waves (continuous red) and the Belle data on the differential 
yy —> n + n~~ cross-section at \Ts = 0.95 GeV. 



_yy — > n + n /K + K data 

We use the previous fitted values of the parameters from yy — > 
n () 7T° for predicting the yy — > 7i + n~ process. The results for the 
total and differential cross-sections are shown in Fig. |6l where 
one can notice that the prediction is not good between 0.5 to 0.9 
GeV if one only retains the S and D-waves in the partial waves 
of the Born term. A similar problem has been encountered by 
the authors of l34ll (Fig 4) who solve this (Fig^R)]) by introducing 
a form factor (which looks ad hoc) from [ 36] for the D-wave 
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Born contribution which is important in the charged channel PI 
Instead, we improve our analysis by adding all higher partial 
waves in the Born term and by including the direct terms of the 
S -waves. One can see in Fig. [6], that the predictions using the 
parameters from yy — > n rn° are quite good compared with the 
data without introducing any ad hoc form factors. 

8. On some possible substructures of the cr and /o(980) 

Our work has been motived for shedding light on the possible 
substructure of the cr// (600) and /o(980) by analyzing their 
yy widths, where the value of direct width is of special interest 
for a comparison with theoretical calculations based on quark 
and/or gluon loops. 

- The cr meson: in our previous analysis of the hadronic cou- 
plings from nn — > nn, KK processes, we have noticed 
that because of the large size of its coupling to nn and KK 
[\g<rK+K-\l\gan+n-\ = 0.37(6)] (see also JSHli), the cr cannot 
be mainly a nn molecule or/and a four-quark state components, 
while its broad width cannot be explained by a large qq compo- 
nent. These observations go in line to a large gluon component 
expected from a low-energy theorem analysis in flo[[l2i fl3l[^l 
The averaged result of the yy direct width of about 1 .2(3) keV 
(see Eq. [63) , when runned at the on-shell <x mass, can indicate 
an eventual large gluon component of the cr Moreover, our 
result seems to rule out a large four-quark component which 
has a too small direct yy width of about some few eV [41] from 
QSSR or from some cr-like models [42]. 

- The /o(980): the direct yy width of 0.27 keV for the / (980) 
in Eq. doTt is neither compatible with a large four-quark com- 
ponent (too small) nor with a large S 2 component (too big). A 
pure is component where a yy width is expected to be about 0.4 
keV Il2l 1 1 3L 14311 is not however favoured because of the non- 
zero coupling of / (980) to nn [\g fK *K-\/\gf^\ = 2.59(1.34)] 
obtained in HH. 



9. Gluonium production from J /iff and radiative decays 

The previous possible gluonium assignement of the cr and 
/o(980) can be tested from the J /if/ and (p radiative decays. This 
can be done following the work of [44], using dispersion rela- 
tion techniques and the Euler-Heisenberg effective Lagrangian, 
where the gluonic part of the amplitude can be converted into 
the non-perturbative matrix element (0\a s G 2 \S ) ~ M 2 s f$ , where 
the decay constant fs (S = cr B , G) can be obtained from QSSR 
which predicts two scalar gluonia {thecrsil) coupled strongly 



15 A similar observation can be done for the process yy — > K + K~, where 
the Born contributions due to the S and D-waves are also very important and 
dominates over the scalar meson contributions, as can be shown in Fig. 10 of 
13711 . which we have checked. 

16 A large gluon component of the cr meson with a Breit-Wigner mass and 
width of about 1 GeV has been also advocated in 1 17]. A large gluon component 
for the <x and /o(980) also emerges from a calculation of the gluonia spectrum 
using a Nambu-Jona-Lasinio model 11811 . 

17 An isoscalar S 2 = ( 1 / V2) (uu + dd^j state is expected to have a yy width of 
about 4 keV GJ.QJ] from QCD spectral sum rules (QSSR) [38, 39] and quark 
model Bch . while a pure glu onium state <xg with M„- B = 1 GeV has a width of 
about (0.2-0.6) keV ll(j[lj from some low-energy theorems. 



to nn, KK due to OZI violations and the G( 1 .5- 1 .6) coupled to 
i]'r]' ', rpf obtained using quenched lattice) |T3, 12, 13]. In this 
way, one obtains: 



r(j/iff ->sy) 



(M. 



^656100 \ M c 
r(J/i//^> e + e-) 



Ms 



fs 



(66) 



where fi\ = -1/2(1 1 - 2n//3) for S f/(n)/ is the first coefficient 
of the B function; M c is the constituent charm quark mass which 
we take to be about Mj/^,/2. Using the decay constant fs = 
f aB a (1.4 ~ 1.0) GeV for M ff8 * (0.75 - 1) GeV ESElEl, 
one predicts: 



B(J/\p -> o- B y) x B(o~b -> all) * (0.4 ~ 1.0) x 10~ J , (67) 



which one can compare with II23I1 : B(J/if/ — > rj'y) (4.7 + 
0.27)xl0~ 3 and B(J/iff -> f 2 y) =* (1.43+0.1 l)xl0~ 3 . Applying 
this result to the case of the G(l .5- 1 .6) glueball of higher mass, 
one obtains II 1 0L 1 1 2L 1 1 311 : 



B(J/if/ -> Gy) (5.0 ± 3.8) x 10 



-4 



(68) 



where we have used fc — (0.39 + 0.15) GeV. This prediction, 
which needs to be improved, is in agreement with the experi- 
mental lower bound of (5.7 + 0.8) x 10~ 4 [23]. Extrapolating 
this analysis to the case of the 0-meson, one obtains by using 
the strange quark constituent mass M s M^/2: 



B{(j)^Sy)^ 1.2x10" 



(69) 



which, despite the crude approximation used, compares quite 
will with the data Q: B((f> -> n + n~y) a (0.41 + 0.1) x 10" 4 
and B(<p /o(980)y) ^ (3.22 + 0.19) x 10" 4 . 



10. Gluonium production from D s semi-leptonic decays 

An analogous analysis has been done for D s semileptonic de- 
cays, where one expects that the gluonium production will be 
of similar strength as the one for a qq state [45]: 



rjA -> <r B (gg)lv] 
F[D S -> S 2 (qq)lv] 



1 



fa- B 



www IS' * 0(l> 



(70) 



for fo- B x 1 GeV, where |/ + (0)| ^ 0.5 JH is the form fac- 
tor associated to the qq semileptonic production. Some other 
hadronic processes like J/tf/ — > a>, <p + nn/KK (see e.g. [5]) and 
D( S ), — > 3n, ... (see e.g. B47I1 ) decays can also be studied, but 
the mechanism for the decays are expected to be more complex 
than the one of radiative and semi-leptonic processes discussed 
previously. However, these analysis emphasize the important 
role of the kaon loops contributions, which in our approach are 
due to the rescattering contributions. We plan to come back to 
these issues in a future work. 
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11. Conclusions 

At first sight, a gg-gluonium mixing scheme like the one pro- 
posed below 1 GeV in [ 12, 48] might be appropriate for describ- 
ing the cr and /o(980) obtained from our fits of nn — > tttt/KK 
and yy — > nn scatterings. The values of their "direct widths" 
favour a large gluon content for the cr meson but are not deci- 
sive for explaining the substructure of the /o(980) meson. How- 
ever, the large values of the rescattering widths, due to meson 
loops because of the large couplings of the <x and /o(980) to 
nn or/and KK, can be also obtained if they are gluonia states 
but not necessarily if they are four-quark (diquark-antidiquark) 
or molecule states as currently claimed in the existing litera- 
ture. The agreement of our predictions for a gluonium produc- 
tion through radiative <p radiative decays with the data seems to 
support some large gluon component for the cr and to a lesser 
extent for the /o(980). This test can be pursued in the analysis 
of J /iff radiative and D s semi-leptonic decays. We plan to ana- 
lyze in details the substructure of these light scalar mesons by 
including mixings in a future work. 
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